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We have performed numerical simulations of the unpolarized e'^e~ — > pp process in kinematic 
conditions under discussion for a possible upgrade of the existing DAFNE facility. By fitting the 
cross section angular distribution with a typical Born expression, we can extract information on the 
■ ratio \Ge/Gm\ of the proton electromagnetic form factors in the timelike region within a 5-10% 

' uncertainty. We have explored also non-Born contributions to the cross section by introducing a 

' further component in the angular fit, which is related to two-photon exchange diagrams. We show 

that these corrections can be identified if larger than 5% of the Born contribution, and if relative 
' phases of the complex form factors do not produce severe cancellations. 
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I. INTRODUCTION 



> 

0^ ' The form factors of hadrons, as obtained in electromagnetic processes, provide fundamental information on their 
internal structure, i.e. on the dynamics of quarks and gluons in the nonperturbative confined regime. A lot of data for 
nucleons have been accumulated in the spacelike region using elastic electron scattering (for a review, see Ref. jjj and 
\^ [ references therein) . While the traditional Rosenbluth separation method suggests the well known scaling of the ratio 
Ge/Gm of the electric to the magnetic Sachs form factor, new measurements on the electron-to-proton polarization 
transfer in e~p — *■ e~p scattering reveal contradicting results, with a monotonically decreasing ratio for increasing 
momentum transfer = —Q^ M- This fact has stimulated a lot of theoretical work in order to test the reliability 



I ' of the Born approximation underlying the Rosenbluth method (see Ref. |^ ^ |^ and references therein) , and made it 
Oh critical to deepen our knowledge of G^; and Gm also in the timelike region by mapping the dependence of their 
^ moduli and phases. 

'T"! , Timelike form factors, as they can be explored in e'^e~ — > HH or pp £~^£~ processes, are complex because of the 
residual interactions of the involved hadrons H (protons p). Their absolute values can be extracted by combining the 
measurement of total cross sections and center-of-mass (cm.) angular distributions of the final products. The phases 
are related to the polarization of the involved hadrons (see, e.g., Refs. BHi). 

Experimental knowledge of the form factors in the timelike region is poor (for a review see, for example, Ref. Q). 
There are no polarization measurements, hence the phases are unknown. The available unpolarized differential cross 
sections were integrated over a wide angular range, such that the relative weight of \Gm\ and \Ge\ is still unknown. 
In the data analysis, either the hypothesis Ge = or \Ge\ = \Gm\ were used. While the former is arbitrary and 
certainly wrong near the physical threshold — 4m^, with m the nucleon mass, the latter is valid only at q^ = Am? 
but unjustified at larger q^ . As for the neutron, only one measurement is available by the FENICE collaboration 
for < 6 GeV^, which displays the same previous drawback. 

Nevertheless, these few data reveal very interesting properties. The amplitudes in the timelike and spacelike regions 
are connected by dispersion relations In particular, \Gm\ should asymptotically become real and scale as in the 
spacelike region. However, a fit to the existing proton \Gm\ data for q^ < 20 GeV^ is compatible with a size twice 
as larger as the spacelike results Moreover, the very recent data from the BaBar collaboration on \Ge/Gm\ [l^ 
show that the ratio is surprisingly larger than 1, contradicting the spacelike results with the polarization transfer 
method 0] and the previous timelike data from LEAR Also the few neutron data for \Gm\ are unexpectedly 
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larger than the proton ones in the corresponding range Finally, all the available data show a steep rise of \Gm\ 
at q^ Am?, suggesting the possibility of interesting (resonant) structures in the unphysical region (for more details, 
see Ref. [H). 

The possible upgrade of the existing DAFNE facility 0, to enlarge the cm. energy range from the (p mass to 
2.5 GeV with a luminosity of at least 10"^^ cm^^s"^, would allow to explore with great precision the production of 
baryon-antibaryon pairs from the nuclcon up to the A. 

In the following, in Sec. ^1 we briefly review the formalism necessary to extract absolute values and phases of baryon 
timelike form factors from cross section data. In Sec. IIIII we outline the main features of our Monte Carlo simulation 
for the measurement of proton form factors. In Sec. IIVI we focus on unpolarized pp production, and in particular on 
the relevant problem of a precise extraction of the ratio = \Ge/Gm\ from the cos^ 9 term of the expected event 
distribution at any given q^. The single-polarized case will be discussed elsewhere Sec.0 we apply the same 

procedure to identify the presence of a contamination in the unpolarized event distribution by processes beyond the 
Born approximation, in particular contributions from two-photon (27) exchange. Some concluding remarks are given 
in the final Sec. EH 



II. GENERAL FORMALISM 



The scattering amplitude for the reaction e+e^ BB, where an electron and a positron with momenta ki and fc2, 
respectively, annihilate into a spin-i baryon and an antibaryon with momenta pi and p2, respectively, is related by 
crossing to the corresponding scattering amplitude for elastic B scattering. In the timelike case, and in presence 
of a possible 27 exchange term, quite a few independent real scalar functions are required to describe the process. 
There are several equivalent representations of them; here, we use the one involving the axial current following the 
scheme of Ref. The scattering amplitude can be fully parametrized in terms of three complex form factors: 

GE{q^,t), GM{<fit), and 0^(9^, t), which are functions of q^ = {ki + ^2)^ and t = (fc2 — PiY ■ They refer to an 
identified baryon-antibaryon pair in the final state, rather than to an identified isospin state. In the following, we will 
use cos 6 instead of the variable t, with 6 the angle between the momenta of the positron and the recoil proton in the 
cm. frame. 

In the Born approximation, Ge and Gm reduce to the usual Sachs form factors and they do not depend on cos0, 
while Ga ^ 0. For the construction of our Monte Carlo event generator, we rely on the general and extensive relations 
of Ref. |l8| up to the single polarization case. These relations assume small non-Born terms and include them up 
to order (with a the fine structure constant), i.e. they consider 27 exchange only via their interference with the 
Born amplitude. These corrections introduce explicitly six new functions, that all can depend on cos0: the real and 
imaginary parts of AGb and AGa/, i.e. the 27 corrections to the Born magnetic and electric form factors, and the real 
and imaginary parts of the axial form factor Ga- The role of these corrections in the unpolarized and single polarized 
cross sections is very simple: it can be deduced from the Born term by adding the contribution of Ga and substituting 
the Born form factors with the "27— improved" form factors, i.e. GE,M{<f) Ge,m{(1^) + ^GE.M{<f ,cos9). In 
Ref. :18j, this fact is somehow hidden by neglecting terms of order a*. However, for sake of simplicity we keep 
27 effects via the axial form factor only. Within this scheme, when summing events with positive and negative 
polarization, the unpolarized cross section can be written as 



= a{q')[l + R{q') cos' e]-b{q^)Yle[GM{q')GA*{q\cose)] cosO , (1) 



dcos9 

= #f - (-IGmP + \GEn , b{q') ^ ^ , (2) 

2q^ T \ T ^ ' q'^ T 

2^ ^ r\GMiq')\'-\GEiq')\' ^ _l_ 

^ r\GM{q'W + \GE{q'W ' 4m2 ' ^' 



Measurements of the unpolarized cross section at fixed q' for different values of 9 allow us to fit the different cos" 
terms, from which we can extract \Ge/Gm\ and some information on GA{q^ ,cos9). 

For spin-i baryons with polarization Sb, the cross section is linear in the spin variables, i.e. dcr — da° (1 
with d(7° from Eq. and A the analyzing power. In the cm. frame, three polarization states are observable 
the longitudinal Vz, the sideways Vx^ and the normal Vy. The first two ones lie in the scattering plane, while the 
normal points in the pi x k2 direction, the x,y,z, forming a right-handed coordinate system with the longitudinal 
z direction along the momentum of the outgoing baryon. The is particularly interesting, since it is the only 
observable that does not require a polarization in the initial state |a, 0| ■ With the above approximations, it can be 
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deduced by the spin asymmetry between events with positive and negative normal polarizations: 

Vy - j^ rf^I"^^! - ^^^0[^o.eiu,[GM{q')GW)\ - ^^Im [G^(g2) G^*(g^ cos0)] } . (4) 

This spin asymmetry can be nonvanishing even without polarized lepton beams, because it is produced by the 
mechanism pi x k2 • Ss, which is forbidden in the Born approximation for the spacelike elastic scattering 't^. However, 
the measurement of Vy alone does not completely determine the phase difference of the complex form factors. By 
defining with 5e and 5m the phases of the electric and magnetic form factors, respectively, the Born contribution 
is proportional to sin(JM — 5e), leaving the ambiguity between [5m — 5e) and tt — {5m — 5e)- Only the further 
measurement of Vx can solve the problem, because Vx oc Re(Gj^^ G^) cx cos((5m — 5e) 0- But at the price of 
requiring a polarized electron beam. 

III. GENERAL FEATURES OF THE NUMERICAL SIMULATION 

We consider the e+e^ pp process. Events are generated in the variables , 9, 4> (the azimuthal angle of the 
proton momentum with respect to the scattering plane), and Sy (the proton polarization normal to the scattering 
plane). Then, the distribution is summed upon the spin and integrated upon 0. 

The exchanged timelike is fixed by the beam energy. Only the scattering angle 9 is randomly distributed. For a 
given q^ , the Born term is responsible for the [1 + R{q'^) cos^ 9] behaviour of the event distribution. Since at this level 
and with the discussed approximations no other dependence in 9 is present, observed systematic deviations from the 
[1 + R{q^) cos^ 9] behaviour will be interpreted as a clear signature of non-Born terms. 

An overall sample of 300 000 events has been considered with 3.8 < q^ < 6.2 GeV^ and |cos0| < 0.9. Since the 
integrated cross section for e+e~ pp in the considered region is approximately 1 nb, at the foreseen luminosity of 
10'^^ cm~^s~^ this sample can be collected in one month with efhciency 1. The upper q^ cutoff is consistent with the 
upper limit of the cm. energy in the presently discussed upgrade of DAFNE The lower limit includes the 

pp threshold, but in our simulation we avoid this region which is characterized by peculiar mechanisms like, e.g., the 
Coulomb focussing and possible subthreshold resonances. 

Results depend much on the binning of the events. In fact, while a larger number of bins in the same {q^,9) 
region implies a more precise determination of q'^ and 9, at the same time it leaves each bin with much less events and, 
consequently, with larger error bars. We find the best compromise with 6 equally spaced q^ bins with width Ag^ = 0.4 
GeV^, and 7 equally spaced cos 6* bins with width Acos0 w 0.257 for the solid angle | cos 6*1 < 0.9. Incidentally, we 
notice that with Ag^ = 0.4 GeV^ the deviation inside each bin between the statistical average q^ , calculated in the 
simulation, and the central value (q^ax ~ Imin)/'^ about 3%. 

Since it is well known that |Gm(q^)| ~ l/'Z'* E3|> from Eqs. (|1I3|I we expect the cross section to approximately 
fall like 1/q^^. This implies that bins at higher are scarcely populated. For example, out of the 300 000 selected 
events only 700 are accumulated in the 5.8 < q^ < 6.2 GeV^ range for each cos 6* bin. In these conditions, a gaussian 
fluctuation of la differs by 5% from the average value; therefore, deviations of at least 10% from the Born cross 
section are needed to be recognized as systematic effects due to non-Born terms. The situation is evidently more 
favourable at lower q'^: in the 3.8 < q^ < 4.2 GeV^ range around 23 000 events are accumulated in each cos6' bin. 
Of course, it could be possible to asymmetrically split the beam time in order to make the statistics of each q^ bin 
more homogeneous. More importantly, since in a e~^e~ collider experiment q^ is fixed by the colliding particle energy, 
it could be more convenient to concentrate on narrow subranges like, e.g., 4 < < 4.1 GeV^ and 4.4 < q^ < 4.5 
GeV^. The reconstruction efficiency does not change much when the same number of events is concentrated in two 
bins centered around the same q^ value but with different widths Ag^, at least for Aq^ <C 1 GeV^. We define as 
"default conditions" the ones where the whole possible range in is divided in equally spaced bins and each one is 
sampled with the same beam time; our results are presented in this framework. A discussion about other possible 
choices is beyond the scope of this paper. 

Events can be generated only by inserting specific parametrizations of the proton form factors in the cross section. 
For Ge and Gj\f , several models can be considered in the timelike region 7, 8], mostly derived from extrapolations from 
the spacelike region. We selected the parametrizations of Refs. ,20, 21], because they have been recently updated 
in Ref. Q by simultaneously fitting both the spacelike and timelike available data. Moreover, both cases release 
separate parametrizations for the real and imaginary parts of Gb and Gm- We indicate the former as the IJLW 
parametrization (from the initials of their authors), and the latter as the Lomon parametrization. 

Since the explored q^ range is not large, it is reasonable to make also much simpler, but equally effective, choices. 
Accordingly, we use also a parametrization where all form factors are real and proportional to the same dipole term 
1/(1 -I- q'^/ql)^, with qo — 0.71 GeV, as in the spacelike region. This choice looks reasonable also for Ga, of which 
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nothing can be said except that hopefully its modulus could have the same dipole trend in the considered range. 
All form factors being proportional to the same function, the actual parameters are the ratios r,, — \Ge/Gm\ and 
Ta = \Ga/Gm\, and the relative phases /3e and /3q of and Ga, respectively, with respect to Gm\ they determine 
the sign of the corresponding ratios. We will refer to the parametrization DipO when = Va = 0; to Dipl when 
Te = 1 and ra — 0; to Dip27 when Vf. — 1 and = 0.2. Since the Born term in the cross section depends only on 
IGbP and |GmP, we can safely take fie — 0. We refer the reader to Sec.^for a discussion about ra and Pa- 
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FIG. 1: Angular distribution of 300 000 events for the e^e~ — » pp process at 3.8 < < 6.2 GeV'^ and | cos6'| < 0.9, and with 
the Dipl parametrization of proton form factors (see text). The points correspond to 7 equally spaced cosO bins for each of 6 
equally spaced bins. Solid curves are the results of a 3 parameters angular fit with B cc ra = \Ga/Gm\ = (see text). The 
highest curve corresponds to the lowest bin 3.8 < < 4.2 GeV'^; next lower curve to the adjacent bin 4.2 < < 4.6 GeV'^, 
and so on. 



IV. RECONSTRUCTION OF \Ge/Gm\ 

As already stressed in the Introduction, one of the puzzling features of the available data is that the BaBar 
collaboration reports a value for \Ge/Gm\ bigger than 1 in the timelike region which contradicts the spacelike 
results and the previous timelike data from LEAR Therefore, it is important to explore the possibility of new 
and more precise measurements of \Ge/Gm\ in the unpolarized e+e~ — > pp process at the upgraded DAFNE [TslIT^. 

Our analysis consists of three steps; for sake of simplicity, we describe it in the following for the case of the 
parametrization Dipl: 

1) We generate 300 000 events for 3.8 < < 6.2 GeV^ and |cos6'| < 0.9 (i.e., for 25° < 6* < 155°). The sorted 
events are divided into 6 equally spaced bins of width Ag^ = 0.4 GeV^. In each q^ bin, the events are divided 
into 7 equally spaced cos6' bins with width A cos 6* w 0.257. In Fig. ^ the points show the 6 distributions in q^, 
each one consisting of 7 points describing the distribution in cos 9. From the discussion in the previous section, 
the highest cos 61 distribution corresponds to the lowest bin 3.8 < < 4.2 GeV^; the next lower distribution to 
the adjacent bin 4.2 < q^ < 4.6 GeV^, and so on. 

2) For each q^ bin, the cos 6 dependence of the event distribution is fitted by a function of the form 

A^(cose')= a(i - B cos9 + R cos^ , (5) 

with three fitting parameters A, B, and R. The results of each fit are represented by the solid curves in Fig. ^ 
The parameter A is not crucial for the following discussion and it will not be considered. The parameter B 
contains the effects of 2j exchange which will be discussed in Sec. [3 Here, we just notice that, with the 
input ra = \Ga/Gm\ = and discarding 2j effects in G^ and Gm, the simulation produces almost perfectly 
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FIG. 2: Fitted R values from the event distribution of Fig.Q according to Eq. Q (see text). Solid line represents the expected 
values for the Dipl parametrization according to Eq. @ (see text). 

symmetric curves in cos 6*, indicating that within the statistical uncertainty _B ~ 0, as expected. The parameter 
R must obviously be identified with the function R{q^) of Eq. (jSjl, and it is related to Ve = \Ge/Gm\ by 



R = 



(6) 



The fitted R values are reported in Fig. |3 against the corresponding average values of each bin. The error 
bars are the output of the fitting procedure, i.e. they represent the uncertainty in determining R when using 
the form to fit the event distribution in Fig.^ The solid curve in Fig. [3 represents the expected value of R 
when the Dipl parametrization for the form factors is inserted in Eq. Q through Tg. 
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FIG. 3: Extracted — \Ge/Gm\ values from the event distribution of Fig.Q Solid line gives the expected value according to 
the Dipl parametrization (see text). 



3) If we define AR as half of the error bar associated with a given R value in Fig. [3 we can also define the 
pair {Rmin = R — A_R, Rmax = R + AR). By inverting the relation it is possible also to deduce the 
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corresponding pair {r^^^^, r^^^^), which can be transformed into {r^, Are) with Ar^ = ('"e„„^ — ^e„i„)/2. These 
extracted values of = \Ge/Gm\ and the associated error bars Ar^ are displayed for each corresponding 
bin in Fig. |31 The solid line represents the expected value (in this case, re = 1 because we are using the 
Dipl parametrization) . The comparison between expected and reconstructed re values displays the potential 
reliability of the measurement of with a sample of 300 000 events and with the chosen binning in . 

Our reconstruction procedure is based on the structure of the cross section which naturally suggests to fit the 
directly observable parameter R and then to link it to according to Eq. 0. It must be remarked that the usual 
formula for the error propagation cannot be used to derive Are from AR, because for some R values the relative error 
is very large (see Fig.l^))- We stress that the error bars are related to the fitting procedure, i.e. they give the integral 
deviation between the test function © and the data points. They are not directly related to the statistical errors on 
the average population of each bin. In the measurement of a really random variable, in the limit of infinite number 
of events, the statistical error on the population of a single bin tends to vanish because of the central limit theorem. 
In the same limit, the error from the fitting procedure does not tend to vanish unless the number of parameters is 
equal to the number of bins (in this case, 7 cos 9 bins for each independent fitting procedure). Anyway, the statistical 
fluctuations can be estimated from the distance between expected and reconstructed values. 

In order to extract r^. from R we need to assign a common r value, i.e. g^, to all the events falling into a specific 
bin. Integrating over all the generated events, we have calculated the average q"^ for each q'^ bin, and used it for the 
extraction of re. We recall that this average differs by about 3% from the central value {q^^x ~ Imin)/'^ deduced 
from the bin width Ag^ = g^^^ _ g^^^. 

It is important to notice also that, despite the large number of available events, the region with the smallest is 
delicate. In fact, R is small but has not a small derivative in q^ near threshold. This reflects in a striking contrast 
between the very small fitting error bar, on one side, and a worse comparison between reconstructed and expected 
than at larger g^, on the other side. 
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FIG. 4: Extracted re — \Ge/Gm\ values from the event distribution in the same conditions as in Fig. Q but for the IJLW 
parametrization of Ref. (2C| . as updated in Ref. Q (see text). Solid line is the expected result. 

We have repeated the reconstruction procedure also for the DipO parametrization, i.e. for = \Ge/Gm\ — 0, which 
has been used in some experiments to extract the timelike \Gm\ data from unpolarized cross section measurements. 
The comparison between expected and reconstructed r^. is more difhcult, because not only the fitting error is obviously 
larger, but also the statistical one. However, since the considered q^ are rather close to the physical threshold, \Ge/Gm\ 
may not be very different from its threshold value of unity. According to the results with the Dipl parametrization, 
this should simplify the task of a precise extraction of in a real experiment. 

To verify this conjecture, we have repeated the procedure starting from the realistic parametrizations IJLW (Fig.^ 
and Lomon (Fig.jSJ. Again, the solid curves show the expected values of re as calculated for each different choice. 

Three peculiar features distinguish these two examples from the previous ones: Ge and Gm are complex objects; r^. 
is not a constant; it is close to 1 but bigger than 1, in agreement with the recent findings of the BaBar experiment |l2j . 
In both cases, the quality of the reconstruction of is very good, confirming that in most of the q^ bins the extraction 
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FIG. 5: Extracted re = \Ge/Gm\ values from the event distribution in the same conditions as in Fig. Q but for the Lomon 
parametrization of Ref. [2H. as updated in Ref. (see text). Solid line is the expected result. 

of \Ge/Gm \ is possible with an accuracy within 10% for the considered sample and binning. This is mainly due to the 
simple dependence of the cos^ term in the angular distribution upon only one parameter, i.e. = \Ge/Gm\', to the 
neglect of more complicated angular terms, like CDs'* 9; and to the vicinity of to the physical threshold = Am? , 
which implies \Ge/Gm\ ^ 1- 

V. THE TWO-PHOTON CONTRIBUTION 

As already anticipated in Sec. ^ we account for non-Born contributions to the unpolarized cross section only 
through the axial form factor Ga- According to the Dip27 parametrization discussed in Sec. II VI Ga is modeled as 

GA{q\cosd) - Tae'^-^ GM{q^) , (7) 

where we neglect any explicit dependence on cosO. We take Gm real and with the dipole form. Indeed, in the 
kinematical region under analysis the ratio Ge/Gm is complex, but most likely dominated by its real part. A 
nonvanishing phase is due to the reinteraction between the two hadrons in the final state. Very close to the threshold 
(i.e., in a relative s wave) the phenomenon presents peculiar features, but from the onset of the p wave contribution 
onwards, an approximately semiclassical black-disk regime takes over 22]. In these conditions, the main effect of the 
pp interaction is a flux damping. The transition from the near-threshold to the black-disk regime is signalled by the 
change of sign of the p parameter, i.e. the ratio between the real and imaginary parts of the forward pp scattering 
amplitude. Such a transition takes place below the q^ range considered in our analysis. In the black-disk regime, it is 
natural to think that the flux absorption leads to the appearance of an imaginary part for Ge and Gm- However, the 
absorption cross section is much below the unitarity limit, suggesting that the phases of Ge and Gm may be small. 
For Ga it is not possible to conclude that it is " quasi- real" . |Gyi| is presumably small, and the cuts on the two-photon 
lines involving all the possible on-shell intermediate states in the amplitude, may produce a relevant imaginary part. 

The Ga form factor enters the unpolarized cross section with a contribution proportional to cos0Re(Gj(jG^). 
Therefore, we can have 

Re[GM G\] « Gm MG*a] ~ Va cos /3,^ G^, = rr, GIj . (8) 

In our simulation with the Dip27 parametrization, we consider Gm, Ge-, and Ga-, all real functions of q^, with 
Te = \Ge/Gm\ = 1 and rr^ k Va — \Ga/Gm\ ~ 0.2, i.e. with 13a — 0. This is the most optimistic situation we may 
guess, taking into account that: (i) no realistic models are available for Ga] (ii) in the spacelike case, the assumption 
that ra is at most 0.2 does not contradict the data; (iii) in the timelike case, Ga should have a nonvanishing phase. 

In Fig. El the generated events are sorted in q^ and cos6' bins as in Fig.^ We can see now that 27 effects show up 
in a marked left-right asymmetry of points in cos 9. The solid curves are obtained from the fitting formula (jSJ, where 
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FIG. 6; Angular distributions for 300 000 events for the e^e pp process in the same conditions as in Fig. Q but for 
ra = \Ga/Gm\ =0.2. 



now the parameter B is given by 



with 



f{q\r,) =4- 



(9) 



(10) 



Fitted B values | 




FIG. 7: Fitted B values from the event distributions of Fig. [^according to Eq. ^ (see text). Solid line represents the expected 
value for the Dip27 parametrization according to Eq. JUJ (see text). 

In Fig. 13 the fitted values of B and the relative errors are shown. As usual, the solid line shows the expected result, 
by inserting — 1 and ra = 0.2 in Eqs. H9I10|I . 
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In order to estimate from B, we assume small errors for both B and r^, with the caveat that such an assumption 
is not valid near the lowest threshold. Using the formula for the systematic error propagation, we get 



Ar - — 



AB 



■BAr, 



(11) 
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FIG. 8: Extracted ra = \Ga/Gm\ values from the event distributions of Fig. |Sl according to the Dip27 parametrization; solid 
line is the expected result (see text). 

In Fig.lHl the Va values, reconstructed from Eq. are shown together with the relative errors Ar^ from Eq. 
The solid line is the expected value ra = 0.2. From the figure it is evident that the error bars (fitting errors) and the 
discrepancy from the expected value (statistical errors) are of magnitude 2Ara ~ 0.05. Therefore, with the considered 
sample and binning we would be able to distinguish a non-Born contribution when it is of the size \ra\ > 0.05. 

The additional non-Born corrections to and Gm can be effectively reabsorbed in the measurement of the 
left-right asymmetry in the cos 9 distribution. As shown in Ref. charge conjugation invariance imposes general 
symmetry properties of the Born and 2j amplitudes with respect to the cos 9 — cos 9 transformation. In particular, 
the 27 corrections to the form factors should respect the following constraints 



AG£;^A/(g^cos6') = -AGi^,M(g^-cos6») 



Ga [q^ , cos 61) = Ga (g^ - cos 9) 



(12) 



It follows that the contribution of the Born < 
expression 



27 interference term to the unpolarized cross section has the general 



dcos9 



cos6'[co(g^) -f ci(g^)cos^6l -I- C2(g^)cos^6l + ...] 



(13) 



where Ci {i = 0,1,...) are real coefficients incorporating effects from all the three different form factors. In our 
analysis, where we take AGe,m = and we consider only the first term in the power expansion in cos9 of Ga, the 
left-right asymmetry proportional to gives us only a lower bound for the actual strength of the 27 effects. It is 
however evident that several independent observables, including the polarization of the recoil proton and/or electron 
beam, are necessary to disentangle the two-photon contribution from each of the three form factors. 



VI. CONCLUSIONS 



We have performed numerical simulations of the e+e ^ pp unpolarized process using a sample of 300 000 events in 
the kincmatical region with 3.8 <q^< 6.2 GeV^, distributed over 6 equally spaced bins. For each q^ bin, events were 
further distributed over 7 equally spaced bins in cos6', with | cos 6*1 < 0.9. For each q^ value, the distribution in cos9 
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was fitted by the function A[l — B cos 9 + Rcos^ 9], where the fit parameters R and B allow for the reconstruction 
of the underlying values of the ratio = \Ge/Gm\ and = \Ga/Gm\, respectively, once models for the proton form 
factors are inserted as input. 

Concerning r^, the reconstruction seems to be relatively simple if G^; ~ Gj\/, as we expect in the considered 
kincmatical range. The size of the sample reproduces the expected \Ge/Gm\ within 5-10%. The worst performance 
is for Ki 6 GeV^, where the bins are more scarcely populated. 

The additional contributions related to two-photon exchange diagrams show up in a left-right asymmetry of the 
angular distribution, driven by the parameter B. In our analysis, we modeled such a contribution in terms of the 
axial form factor Ga, providing a lower bound for the actual strength of two-photon effects. The simulation shows 
that it is possible to identify and estimate this term, if \Ga\ (or equivalently lAG^I or \AGm\) is larger than 5% 
of |Gm|, and the relative phases of the form factors do not combine in such a way to produce severe cancellations. 
However, a more refined analysis, simultaneously including also polarization observables, is needed to better constrain 
and disentangle the two-photon contribution from each of the three form factors. 
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